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Galactic dark matter is modelled by a scalar field in order to effectively modify Kepler's 
law without changing standard Newtonian gravity. In particular, a solvable toy model 
with a self-interaction U (<£>) borrowed from non-topological solitons produces already 
qualitatively correct rotation curves and scaling relations. Although rclativistic effects 
in the halo are very small, we indicate corrections arising from the general relativistic for- 
mulation. Thereby, we can also probe the weak gravitational lensing of our soliton type 
halo. For cold scalar fields, it corresponds to a gravitationally confined Boson-Einstein 
condensate, but of galactic dimensions. 



1. Introduction 

\i\t\\ _ 

The existence of dark matter was postulated for the first time by Zwickj in 1933. 
He showed that the virial masses of clusters of galaxies exceed the mass of luminous 
matter by roughly one order of magnitude. In the 1970's spectrographs became 
sensitive enough that rotation curv es of sp iral galaxies could be observed in large 
numbers by Rubin and collaborator J I I They all showed a systematic tendency 
to stay flat as function of galactocentric radius, v v (r) ~ const., whereas the mass 
distributions of the bulges and the radially exponentially declining surface densities 
of the galactic disks lead to the prediction of falling rotation curves. Evidence for an 
extra mass component became even more decisive when Hi-rotation curved^ were 
observed, because the neutral hydrogen gas ext end s radially outwards many disk 
scale lengths. Even though there is a long debate^! to what degree the inner parts 
of galaxies are dominated by dark matter, the flat outer Hi-rotation curves prove 
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beyond doubt the existence of dark matter halos. The analysis of the kinematics of 
satellite galaxieJ^SI showed that the dark halos extend out to several hundred kpc 
with the cumulative mass still rising linearly with galactocentric radius, M(r) oc r, 
which corresponds to flat rotation curves. 

Persic et a/!^ investigated statistical correlations between structural parameters 
of visible and dark matter using a data set of almost 1000 galaxies and developed 
the model of the so called 'universal rotation curve'. Accordingly galaxies of low 
luminosity are more dominated by dark matter than bright galaxies. Dwarf irregular 
galaxies and low surf ace br ightness galaxies seem to be dominated completely by 
dark matter. BurkerlP2Di2l finds for such galaxies an empirical fit to the rotation 
curves which implies asymptotically tc — > lnr/r, i.e. a logarithmic modification of 
the Kepler law. 

For several classes of gravitational theories, it has been argued that the introduc- 
tion of dark matter is necessary 221 , although in the modified Newtonian dynamics 
(MOND) scheme, exact expressions were derived for the Newtonian potential and 
for the rotation curv However, recent observational constraints^ make it more 
and m ore challenging to interpret galactic rotation in terms of new gravitational 
physical 

An i ntriguing possibility is a scalar field model of the halo first proposed by one 
of u J®7I68 I (?9 I 7() I 71 | using a complex massless field coupled to the Einstein equation, 
an idea which was later taken up by several other authors. Hence, the complex 
massless scalar field forming star-like galactic halos embodies another form of dark 
energy. 

Massive compact halo objects, so- calle d MACH OiP, consisting of baryonic mat- 
ter or, alternatively, of boson stars (BSs}^Z1 5Q I 7$ l d n ot seem to be sufficient to 
resolve this problem completely. More recent approacheJ^^ propose connections 
between scalar field matter and the cosmological constant. 

On the other hand, the standard cold dark matter (CDM) model, where dark 
matter particles interact only through gravity, may be in conflict with observational 
featureson Mpc scale. In particular, the CDM model of Navarro, Frenk and White 
(NF Wr^p tends to predict density profiles of d ark halos which are cuspy at the 
centeil^2l Therefore, Spergel and SteinhardJS^ proposed that dark matter is self- 
interacting, an idea which found ramification in Refs. I fill 1321 The evolution of the 
self-interacting dark matter has been recently analy zed in more detaiPO. 

In the same vein, we consider here a mode@^^ where dark matter is modelled 
by a primordial scalar field with a self- interaction J7(<I>). Since the observed ro- 
tation velocities are roughly bounded 3, by v v < 10 _3 c, i.e. are non-relativistic, a 
Newtonian type approximation sufficies. 

The formation of dark halos as a gigantic Bose-Einstein condensation in early 



a The maximum velocity within a rotation curve of almost 400 km/s is found in a Sa galaxjlS^J 
the velocity of 367 km/s at the farest measured point belongs also to a Sa galaxjEH 
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Universe is discussed in Sec. |21 In particular, a boson star like toy model with a 
<I> 6 repulsive self-interaction is knowrP^l to allow in the limiting case of flat space- 
time an exact spherically symmetric soliton solution of the corresponding nonlinear 
Klein-Gordon equation, as is shown in Sec. 01 Simulating the halo by such a non- 
topological soliton (NTS) for the positive range of the potential, yields a Newtonian 
mass distribution which provides a qualitatively rather good fit to the rotation 
curve data of dwarf irregular and low surface brightness galaxies, see Sec. 0] and 
the comparison with observations in Sec.|SJ In Sec. El we indicate the extension to 
the general relativistic formulation and speculated on the superstring nature of the 
dark matter scalar. In Sec. El we display the gravitational lensing properties of the 
NTS halo. 

2. Dark halo as a self-gravitating Bose-Einstein condensate? 

Since Einstein and Bose it is well-known that identical integer spin particles can 
occupy the same ground state. Such a Bose-Einstein condensate (BEC) has been 
experimentally realized in 1995 for cold atoms of even number of electrons, protons, 
and neutrons, see Anglin and KetterkJ^ for a review. In the mean-field ansatz, the 
interaction of the atoms in a dilute gas is approximated by the effective potential 



This leads to a nonlinear Schrodinger equation for the macroscopic matter wave 
field 'J, in this context known as Gross-Pitaevskii equation. In a microscopic ap- 
proach, one introduces bosonic creation and annihilation operators and b, re- 
spectively, satisfying [&,&'] = 1, and finds that every number conserving normal 
ordered correlation function (&i ■ • ■ b n ) splits into the sum of all possible products 
of contractions (bibj) as in Wick's theorem of quantum field theory (QFT). For 
n = 1 one recovers the Gross-Pita evsk ii equation, whereas the next order leads to 
Hartree-Fock-Bogoliubov equationiPSI 

In the context of astrophysics, an assembly of massive real scalar particles has 
already been considered by Bonazzola and PacinP^ in 1966. Later on, solutions 
for "systems of self-gravitating" bosons were considered by Rufimi and BonazzolcP^ 
similarly to the Hartree-Fock atom. Because a real scalar has no anti-particle states, 
the corresponding Klein-Gordon field $ can only be decomposed into positive and 
negative frequency field operators: 



umu = 4i*i 4 - 



(i) 



(2) 



For a spherically symmetric configuration 



— iljJnlt 



nla 




(3) 



nla 
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can be expanded in terms of the wave function of a bound state similarly as for the 
hydrogen atom. Here b n ia and &L;/ a / are again creation and annihilation operators, 
Rf(r) are radial distributions, Y^\e, ip) = (1 /V4n)Pj a] (cos 9) e' 1 ^ the spherical 
harmonics given in terms of the normalized Legendre polynomials, and \a\ < I are 
the quantum numbers of azimuthal and angular momentum. (The nodeless radial 
field solution Rg(r) will later be denoted by P(r).) 
For the ground state of a cold configuration, 

JV 

\N) = \N,0,0,...) :=n & UlO) (4) 

l 

is choseri^, where |0) is a vacuum state in the curved spacetime 'background'. 
The energy-momentum tensor T /JI/ ($) of the scalar field becomes now an operator. 
When coupled to gravity, the vacuum expectation value (T M „) := (iV| : T^ v : \N) 
for the ground state is inserted into the right-hand side of the Einstein equation 
(|24J) . where : : denotes normal ordering of the operator products. 

An excited state \N, n, /, a) := <I> + |0) of positive energy can be obtained by apply- 
ing creation operators. Such a 'gravitational atomic represents a coherent quantum 
state, which nevertheless can have macroscopic size and large mass. The gravita- 
tional field is self-generated via the mean value (T^) of the energy- momentum 
tensor, but remains completely classical, whereas the real scalar is treated to some 
extent as operator. 

Below some critical temperature (jSJ), such configurations of coherent bosons 
would f orm Bose-Einstein condensates on an astrophysical scale^, and following 
T.D. Led^il are called boson stars (BSs). 

Real scalars have the disadvantage that no local symmetry provides a conserved 
particle number, and so one needs to introduce by hand a normalization condition 
in order to stabilize the system. 

It is gratifying to note that BSs with repulsive se lf-in teraction f/(|3>|) considered 
already by Mielke and Scherzei^J and Colpi et o/P^l have their counterparts in 
the effective potential of BEC. Thus a cold BS can be considered as a self- 
gravitating BEC on an astrophysical Moreover, BSs of rather different 
sizes can occur: it could be just a 'gravitational atom'; it could be as massive as the 
presumed black hole in the central part of a galaxy; or it could be an alternative 
explanation for parts of the dark matter in the halo of galaxies, as we are going 
to explain below. BSs, if they exist, would be an astrophysical realization with a 
self-generated gravitational confinement. 

The idea of a cold mini-boson star as a BEC condensate on a galactic scale was 
reiterated, without referring to the earlier papers, in Refs. 021 1251 m a Newtonian 
approximation. A comparison of the Jeans scale of a dark matter halo and the de 
Broglie wave length of the gravitationally confined bosons provides the estimate 
m$ ~ 10 -22 eV/c 2 for a tiny 'bare' mass of the bosons. 

Bosonic particles such as fundamental scalars allow the possibility of a self- 
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gravitating BEC on an astrophysical scale, commonly referred to as BS. However, a 
repulsive self-interaction has the physical effect of delaying the progressive collapse 
of a cosmic assembly formed by accretion. This, in particular, applies to a cold BS 
below the critical temperature 



2nh 2 
mk 



N 



(5) 



where £(z) is the Riemann zeta function, provided an ideal gas approximation is 
employed. However, for the possible formation of a BEC as a dark halo, one has 
to take into account its immersion into the heat bath of the cosmic background 
radiation of nowadays 2.726 K. 

The physical nature of gravitationally coupled scalar fields^l is at best specula- 
tive: From the field-theoretical point of view, a dilaton field tp coupled universally 
to the trace of the energy-momentum tensor via 



is a rather natural candidate. In order to apply this to galactic halos, one needs an 
ultralight scalar field with a Compton wave length A = h/2m of the order of the 
galactic core r c ~ 10 20 m, corresponding to m v ~ 10 -22 eV/c 2 . On the other hand, 
ultra-light 'gravi-scalars' arising in the framework of five-dimensional braneworlds 
are constrained^ by the precisely observed slow-down of the Hulse- Taylor binary 
pulsar due the gravi tati onal wave emission. The universal coupling JBJ provides 
a natural mechanisrrP21 to protect such a tiny mass from renormalization due to 
quantum loops arising in the Standard Model. In order to suppress large long-range 
forces, an approximate global symmetry needs to be postulated which, nevertheless, 
allows a coupling of ip to the Pontryagin term of electromagnetism via (ip/M*)FAF, 
cf. CarrolPl. 

A related proposal is that of WettericlP^ which endows a massless dilaton (as 
a pseudo-Goldstone boson) with a tiny mass via the conformal anomaly. Such 
a dyna mica l scalar "quintessence" may also account for dark energy^. Another 
approach^! to unify clustered dark matter and dark energy is based on a Born- 
Infeld type higher derivative Lagrangian for the same hypothetical scalar field. 

Recently-"*', we analyzed the bifurcations in the mapping of a self-interacting 
scalar field $ minimally coupled to Einstein's general relativity (GR) to a non-linear 
curvature scalar Lagrangian. Intriguingly, the higher-order Lagrangian bifurcates in 
almost linear Einsteinian patches, but with different effective gravitational strength 
and cosmological constant (dark energy) depending on the cosmological scale. More 
precise constraints on dark energy are expected to come from future data for the 
cosmic microwave background (CMB) radiation®, which may also get modified^ 
via the Sunyaev-Zeldovich effect during its long journey. 
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3. Newtonian non-topological solitons 

As a solvable toy model with self-interaction, let us consider the Klein-Gordon 
equation 1251) with a $ 6 type potential 

C/(|$|)=m 2 |$| 2 (l- x |$| 4 ), xl^l 4 < 1 , (7) 

where m is the 'bare' mass of the boson and x a coupling constant, which are 
thought of as constants of nature. The self-interaction in the radial Klein-Gordon 
equation takes the form dU(P)/dP 2 — m 2 — 2>m 2 xP A - In flat spacetime, such a 
model was first considered by one of ua^ for constructing non-topological soliton 
(NTS) solutions. 

For a spherically symmetric configuration and the choice u> = m, the corre- 
sponding nonlinear Klein-Gordon equation simplifies to an Emden type equation 

P" + -P' + 3 X P 5 = , (8) 
x 

familiar from the astrophysics of gaseous spheres. It has the completely regular 
exact solution 



where we introduced the dimensionless radial coordinate x :— mr, and A = 
^/x~P 2 (0) m terms of the central value. The solution depends essentially on the 
nonlinear coupling parameter x, since the limit x would be singular. This 
feature is rath er characteristic for soliton solutions. Already 1978, it has been 
generalized^! to a NTS with (quantized) angular momentum I. 

In the following, it suffices to restrict ourselves to the above given range for 
which the potential t/(|$|) remains positive. If effects of self- gravitation are taken 
into account, this scalar pote ntial ngedsnotto b e bounded from below: In the 
case of general relativistic BS J 44 l ^ l 4 W7 2 l 4 Q I 7Q | 4fl with & se lf_interaction A|$| 4 , it 
has been proven numerically^! that NTS exist even for negative values of A. The 
extens ion to soliton stars b have been considered by Lee and Pang 4 ! as well as 
GleiseJ22 

The canonical energy-momentum tensor of a relativistic spherically symmetric 
scalar field is diagonal, i.e. T M 1/ ($) = diag (p, — p r , —p±, —p±) with 

1 

2 

Pi = P - U , 

P± = Pr - P' 2 , (10) 



P=Uuj 2 P 2 + P' 2 + U) , 



b Our potential J7J is not bounded from below, however, its positive range is sufficient for applying 
our toy model. An example of a bounded potential is f/(|$|) = m 2 |<I>| 2 (l — | < l?| 2 /| ( I > o| 2 ) 2 typical 
for soliton stars, but then simple analytical expressions arc not available and one needs, for large 
radius, to deal with asymptotic solutions like P(r) — > ±a/ Aj (1 + B 2 e 2x ). 
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where ' = d/ dr. The form (|10|l is familiar from perfect fluids, except that the 
radial and tangential pressure generated by the scalar field are in general different, 
i.e. p r =/= p±. This anisotropy of scalar matter should not been ignored, since it 
holds even in flat spacetime, or in the Newtonian approximation. 
From Ijl 0(1 we find in flat spacetime the energy-density 

p = ^- [2P 2 + P' 2 - X P e ] 



Am 2 



2y/x(l + Ah 



A 4 x 2 - A 2 



(l + A 2 x 2 ) 2 



(11) 



(If we would consider a real scalar field instead, there is merely the change of the first 
proportionality constant 2 in the bracket to 1.) The requirement of the positivity of 
[/(<&) at the origin yields the constraint A < 1, which we will adopt in the following. 

The leading term of the Newtonian ty pe m ass concentration 1(11(1 is exactly the 
density law of the quasi-isothermal sphere^ 

P(r) * • (12) 



At large radii the density falls of like p oc r -2 which corresponds to an asymptotically 
fiat rotation curve. Comparing with the quasi-isothermal sphere 1(12(1 . the central 
density of the NTS model is given by po — Am 2 / ' yfx and the core radius is r c ~ 
1/mA. This implies a scaling law for the dark halos of the form 

Po - oc — , (13) 

where A, which may vary from halo to halo, cancels out. 

This nonlinearly coupled scalar field excerts the following radial and tangential 
pressures: 



P 



m 



P' 2 } 



[ X P 6 

P± = ^r[xP 6 -P' 2 ] = 



A 3 7 



Ah 



2VX 



P + Pr 



rn 
~2 

m 2 [P 2 +P' 2 } 
Am 2 



2y/x(l + A 2 x 2 ) 2 
A 3 m 2 (l - A 2 x 2 ) 
2y/x(l + A 2 x 2 f " 2y/x 



A 3 i 



1 



(H 

Thus, at the center of the NTS, we have p r (0) 
radial infinity 



Pr 



-P± 



2y/xAa 



A 2 x 2 ) 2 _ ' " 4l 
p±(0). Asymptotically, we find at 

(15) 



From the energy density ((1111 the mass function 



M{r) 



PV dy 



(16) 
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can be obtained by straightforward integration. With the aid of Reduce!^ we find 
for our model: 



M(r) 



1 



mA^/x 
A „ 



-_m 



6VX 
r 



A 2 -8 
8A 



arctan(Ac) 



A 2 x l + 3A 2 x 2 
~ (1 + A 2 x 2 ) 2 



(2 - A 2 ) 



(17) 



4. Rotation curves 

The tangential velocity v v of stars moving like "test particles" around the center 
of a galaxy is not directly measurable, but can be inferred from the redshift 
observed at spatial infinity, for which 



(1 + Zoo) 5 



1 



(18) 



holds. Due to their non-relativistic velocities in galaxies bounded approximately by 
Um/c < 10~ 3 , we observe z^, ~ v v (as first part of a geometric series) with the conse- 
quence that the lapse function necessarily tends to unity, i.e. e v ~ e Vao — u 2 ) — » 1. 
In general, there are two dominating influences on the redshift of stars: The gravita- 
tional potential of all matter components within the galaxy (stars, gas, dark matter) 
and the Doppler effect. (Since for dwarf galaxies dark matter dominates up to the 
center, we can neglect the contribution of stars and gas.) If e v °° = 1, the redshift 
for receding stars is z = e^ 2 ^ -„,)/(! + v v ) - 1. For example^, the influence 
of the gravitational field of the dark matter on the redshift is just about 1.5 km/s, 
hence, negligible. Of course, if the gravitational potential has the same order of 
magnitude as the Doppler effect, then a velocity of 300 km/s which observations 
reveal so far, would be actually just 150 km/s, but this would certainly demand too 
strong gravitational potentials for galaxies. 

In general, for the static spherically symmetric metric l|27|) considered lateron, 
an observer at rest at the equator of the Schwarzschild type coordinate system 
measures the following tangential velocity squared as a point particle (or a star, 
rega rded as a sufficient small body^l) flies past him in its circular orbit, cf. Misner 
et azE2, p. 657, Eq. (25.20): 



-v 2 i dl P 



-rv' = - [1 - e- x + K Pl r 2 ] e x 



M(r) 2 



r 

M{r) 



exp 



■p T r 



r 



r — kM(t) 



(19) 
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Outside of matter (v = — A) for a weak gravitational field, Eq. I|19|) reduces to the 
Newtonian form t> 2 Nowt ~ M(r)/r. As is well-know n^ ! ? 4 ! a na j ve application of 
the Newtonian limit would have led us to geodesies in flat spacetime, i.e. as if gravity 
would not affect the motion of test bodies like our stars moving in the dark matter 
halo. Thus it is mandatory to go beyond, as is indicated by our approximation of 
the generally relativistic formula l|19|l . Then, also the pressure component p r ^ 
of an anisotropic 'fluid' contributes, our prime example being the case of scalar 
fields. In our NTS model, however, due to the fast decrease of p r , cf. Eq. (|15p. its 
contribution to the asymptotic value of the rotation velocity is almost negligible. 

From the mass function l|17|) of the Newtonian NTS solution and its radial 
pressure (|14|) , we find for the rotation velocity 

9 . 9 (A 2 \ arctanMa;) A 2 A 2 x 2 - 1 , , 

v >~ = 1 + b~ " 1 J + T (TT^W ' (20) 

for which the following approximations 

9 K,A(1 + A 2 )m 2 9 

■* ~ * - im " s 10 " <21) 

hold near the center and at the far field, respectively. Together with l|17fl we con- 
clude that asymptotically M — > 2w^ r. 

Observationally, there is the rough restriction v v /c < 10 -3 of the rotation ve- 
locities of galaxies, which can be used to constrain the mass m < 10 -22 eV/c 2 and 
the coupling constant x °f our NTS model. 



5. Astronomical tests 

The theoretical model predictions against rotation curve data of a set of low surface 
brightness galaxies taken from Ref. |^ El have recently been tested by Fuchs and 
MielkfEHl. There the rotation curves of in total 54 galaxies have been measured 
with high resolution. For about half of them surface photometry is available. For 
these galaxies not only kinematical data have been provided, but also constructed 
dynamical models of the galaxies. The observed rotation curves are modelled as 

v^ir) = i4,buigcM + w 2 !disc (r) + vl Mgas (r) + w 2 !halo (r), (22) 

where v Vt bulge, ?V,disc, f^isgasj and v v ^edo denote the contributions due to the bulge, 
the stellar disc, the interstellar gas, and the dark halo, respectively. The radial 
variations of i> c ,buige(7), u c,disc(0, an d w c,isgas(0 were derived from the observations, 
while the normalizations by the mass-to-light ratios were left as free parameters of 
the fits of the mass models to the data. 

Fits of the form 1)22(1 to observed rotation are notoriously ambiguous. Thus, for 
each galaxy several models are provided, one with zero bulge and disc mass, one 
model with a 'reasonable' mass-to-light ratio of the bulge and the disc, and finally 
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a 'maximum-disc' model with bulge and disc masses at the maximum allowed by 
the data. Furthermore, for each galaxy two types of dark halo models are tried. 
One is the cusped NFW density law and the second is the quasi-isothermal sphere, 
whose observational consequences are recently been reviewed by Bosm^Sl. While 
varying the disc contribution to the observed rotation curve leads to fits of the same 
quality, it is found that the quasi-isothermal sphere models of the dark halos give 
significantly better fits to the data than the cusped NFW density law. Thus the 
scalar field model presented here is in this aspect even superior to the cold dark 
matter model i n it s present form. 

Prada et alW% used data from the Sloan digital sky survey on satellite galaxies 
of isolated host galaxies to probe on 100 kpc scale the outer halo mass distributions. 
They find that the line-of-sight velocity dispersions of the satellites follow closely 
the radially declining velocity dispersion profile of halo particles in a NFW halo. 
This implies an outer mass density distribution of the form p oc r~ 3 which is at 
variance with the prediction of the quasi-isothermal sphere (I12fl . In the cold dark 
matter model the system of satellite galaxi es is assembled during the same accretion 
processes as the dark halo, and Prada et a/!^ assume consistently for the satellites 
the same distribution function in phase space as for the halo particles. In the NTS 
model, however, the dark halo provides for the baryonic matter simply a Newtonian 
force field, for which self-gravity has been neglected in zeroth approximation. The 
distribution function of satellite galaxies in phase space is thus not specified and 
can be modelled according to the observations, even if the potential trough of the 
quasi-isothermal sphere has a shallower profile than a NFW halo. 

Next we have examined the predicted scaling relation A relation of this 

type was found empirically by Salucci and Burkery ^H although this was based 
on the universal rotation curve model of P ersic et a/™2 and the, also empirically 
derived, halo density profile of Burkerl^l This density law resembles the quasi- 
isothermal sphere in that it has also a homogeneous core. In Fig. we show central 
densities versus the inverses of the core radii of quasi-isothermal dark halo models 
constructed in Refs. ^2 El assuming for the discs mass-to-light ratios consistent 
with current population synthesis models. Despite some scatter, there is a clear 
correlation between po and r" 1 over several orders of magnitude, \og(po) oc (1.46 ± 
0.55) log(r f T 1 ). Thus the dark halo model data seem to confirm statistically the 
scaling relation (|13fl . The scatter in the correlation diagram is probably due to 
the near degeneracy of the fits to the observed rotation curves, even if the disc 
contributions to the rotation curves are fixed. Moreover, practice shows that the 
radial exponential scale lengths of the discs cannot be determined more precisely 
than about 20%. This has a considerable effect for the disc contributions to the 
rotation curves IL'l'l) and, consequently, for the dark halo models and might also 
explain some of the scatter of the correlation diagram in Fig.^ We believe, however, 
that this has not changed the general trend. 

Using arguments of the density wave theory of galactic spiral structure, one of 
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Figure 1. Central densities versus the inverses of the core radii of dark halos of low surface 
brightness galaxies derived from modelling their rotation curves. The halo models are constructed 
assuming 'realistic' mass— to-light ratios for the discs. The solid line is the predicted po <x fc 1 
relation. 



us has pointed out, judging from the impli ed intern al dynamics of the galactic discs, 
that the discs might be near to maximurrP^ESESI This would imply for the discs 
of some galaxies mass-to-light ratios which are significantly higher than in current 
population synthesis models of the discs of low surface brightness galaxies. These 
can be modified, though, to yield higher mass-to-light ratios. Therefore we show 
in Fig. |3 the central densities versus the inverses of the core radii of the dark halo 
models, if they are constructed assuming 'maximum discs'. Although the dark halo 
parameters are shifted to other values, the linear correlation persists. With 

log(po) oc (1.08 ± 0.39) log(r- x ) (23) 

the dark halo models fit nearly ideally to the predicted scaling relation (|13|l . In- 
cluded into F ig. HI are also 'maximum disc' dark halo parameters of nearby bright 
galaxieJ^EHl which fit also well to the scaling law. We conclude from this discussion 
that the astronomical tests are rather encouraging for the scalar field model of dark 
halos presented here. 



6. General-relativistic corrections 

In the first part we were focused on a < I >6 toy model. In order to advance to more 
realistic models which include self-gravity, a full general-relativistic framework will 
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be employed which departs from the coupled Einstein-Klein-Gordon equations 

, (24) 
. (25) 

Here R is the curvature scalar, k = 8ttG, the gravitation constant (K = c = 1), 
g the determinant of the metric g^, and $ a complex scalar field. Moreover, 
^ := fVVTsM ®n [yAsls^duJ denotes the generally covariant d'Alembertian. 
The stationarity ansatz 

$(r,i) = P(r)e- lut (26) 

describes a spherically symmetric bound state of the scalar field with frequency 
CO, Note that the case of a real scalar field can readily be accommodated in our 
formalism by putting u) = in our Ansatz. 

For spherical symmetry of the halo, the line-element reads 

ds 2 = e^dt 2 - e x{r Ur 2 - r 2 (d9 2 + sin 2 9d<p 2 ) , (27) 

in which the functions v = v(r) and A = A(r) depend on the Schwarzschild type 
radial coordinate r. 



( dU \ X. 



3, 2008 9:51 WSPC/Trim Size: 9.75in x 6.5in for Proceedings 



main 



13 

Then, the diagonal components of the energy-momentum tensor T^ u ($) — 
diag (p, ~p r , —p±, —p±) generalize to 

p= X - (u 2 P 2 e-» + P' 2 e- X + U) , 
p T = p-U , 

p ± = Pr - P l2 e- X , (28) 
and the Emden equation © to the radial Klein-Gordon equation 

P"{r) + {t^L + H) P '(r) + e x (e" V - P(r) = . (29) 

The decisive non-vanishing components of the Einstein equation are the 'radial' 
equations 

i/ + A' = K{p + Pr )re x , (30) 

A' = npre x - -(e x - 1) . (31) 
r 

Two further components are identically fulfilled because of the contracted Bianchi 
identity V^X^ = which is equivalent to the equation 

^-p T = -v' (p + p T - -(p r -p±)) (32) 
ar \ r / 

of 'hydrostatic' equilibrium for an anisotropic fluid, a generalizatiorP^I of the 
Tolman-Oppenheimer-Volkoff equation. 
The general solution of Eq. I|31|l is 

e- x = l-n^^l-2vl, (33) 
r 

where the mass function M(r) is according to only determined by the energy 
density p. 

According to Eq. (|17|l of the NTS model, the radial metric component e~ x 
asymptotically approaches the value 1 — 2v 2 a < 1. This is not unproblematic: After 
a redefinition of the radial coordinate r — > f := r/^Jl — 2v 2 a , the asymptotic 
space has a deficit solid angle. The area of a sphere of radius r is not Airr 2 , but 
47r(l — 2v 2 xi )r 2 ] cf. analogous results for global monopoles and global texture d*™' 
In order to avoid a linear increase of the mass function which would pose problems 
for the asymptotics, GR allows to redefine the radial coordinate. In the case of the 
more realistic phenomenological Burkert frP^ the velocity tends to zero at spatial 
infinity, with the consequence that no such deficit angle is to be expected. 

So far we have considered a solvable model in flat spacetime. However, when 
the tangential velocity v v = v ip (x) is empirically known, we readily find from (|19(l 
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that the metric components are in general given by 



e" = exp 1 2 J v^dlnxj x zv °° + K, (34) 

e " = iTS" 1 + 2 ' J - (1 -*" 1 ' (35) 

where if is a constant of integration; in the last line we used (|33|) . The approx- 
imation for e x is valid only for p r r 2 ~ 0, i.e. for non-singular radial pressure at 
the origin or for sufficiently fast decreasing pressure at infinity. Fortunately, in our 
NTS model, this condition is satisfied at the origin as well as at infinity, cf. the 
asymptotic function Q15|l. 

For our NTS rotation curve l|20l) . we have derived the scalar field by setting 
the metric components to one in zeroth order Newtonian approximation. Then, we 
calculated the mass function and, consequently, the rotation curve. Now, we can 
do better and determine the metric potentials in first order approximation. For the 
lapse function we find 



exp< 



ln{a; 7 VT 



f 



A 2 x 2 



7/2 



7 'A 

— arctan(Ar) /■ . (36) 



The shift function e A is influenced by both, the rotation velocity 1)21) [I and the radial 
scalar field pressure (|f 4|) . Therefore, a full account of self-gravity seems to demand 
numerical methods. 

A more ambitious approach is to reconstruct a viable scalar potential C/(|$|) 
from the empirical rotat ion curves on the basis of the Einstein equations, similarly 
as in the case of inflatiorM 

As a starting point, we may use the empirical Burkert pronljni which can be 
expressed in the following close analytical form: 



v 



^ B /^ = ^{ln[(f + x) 2 (f + , 2 )]-2arctan(,)} ~ 1 - 



-> 2 . (37) 

x 

ft has a maximum at x — 3.3 in dimensionless units x — r/r c and amounts at 
spatial infinity to a logarithmic modification of the Kepler law. Then, we find for 
the lapse function via Reduce 



e = exp 



1 + x 



2 



x In ' " + 2(1 + x) arctan(x) - In ( 1 + 2x + 2x 2 + 2x 3 + x 4 



(l + x) 2 

(38) 

Assuming a negligible radial pressure term, i.e. p r r 2 ~ 0, the shift function e A = 
1 + 2vfp is simply determined by the Burkett rotation velocity 1)37(1 . 
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For this Burkert fit, we have attempted to reconstruct' a self- interaction po- 
tential U (<&) — p ~ p r from the astronomical data, by resolving the radial Einstein 
equations (|31f) for p and p T with the general result 



U = p-p r 



1 



nre' 



A' 



Accordingly, U is positive when 



d In v v /din x + v 2 + 1/2 > 



(39) 



(40) 



holds, which is satisfied by the Burkert fit. However, since the radial dependence 
of the scalar field $(x) is still unknown, a numerical approach is needed in order to 
reconstruct the potential U($) explicitly. 

In view of the occurrence of scalar fields, like the dilaton or Kalb-Ramond axion 
in effective superstring models leading to axidilaton star\@^, the concordance of 
theory and observations in a viable self-interaction potential deserves further study. 



7. Gravitational lensing of dark halos 

Another promising approach to study the dark halos of galaxies is weak gravitational 
lensing, since it can be measured out to large projected distances from the lense. 
Hence, gravitational lensing of a transparent spherically symmetric NTS dark matter 
halo will be compared here with the scenario of a mini-BS le mP^-U We assume that 
the lens interior is almost empty of baryonic matter, such that deflected photons can 
travel freely through the space between the halo stars which embody only "small" 
disturbanceJ^SI. 

The deflection angle is then given by 

oo 

f be x / 2 

"M = 2 / „ dr - it , (41) 

J r\/r z e v — b z 

ro 

where b — ro exp(— v(ro)) is the impact parameter and ro denotes the closest dis- 
tance between a light ray and the center of our halo. The lens equation for small 
deflection angles can be expressed as 

sin (0-/3) = ^sina, (42) 

-^os 

where D\ s and D os are the distances from the lens (deflector) to the source and from 
the observer to the source, respectively. The true angular position of the source is 
denoted by /?, whereas 9 stands for the image positions. One usually defines the 
reduced deflection angle to be a = 6 — f3 — sin~ x (D\ s sma/D os ) . However, equation 
(|42Jl relies on substitution of the distance from the source to the point of minimal 
approach by the distance from the lens to the source D\ s . For large deflection angles 
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Figure 3. Reduced deflection angle for the NTS dark matter halo. 



the distance D\ s cannot be considered a constant but it is a function of the deflection 
angle so that the form of the lens equation changes^ ^ into 



— — cos cos 

As 



— — sin(# — a) 



[tan0 + tan(a-0)] . (43) 



The reduced deflection angle can be kept defined as a = — [3; cf. Ref. with 
more details on the lens equation. 

Since our dark matter halo, the NTS soliton, has only a weak gravitational field, 
we can apply (|42() together with (|41() . The deflection angle is determined by the 
gravitational potentials (|34|l . (|35|l which are functions of the NTS rotation curve 
(|20|l and the radial pressure l|14|l . As an imput, we use the best fit valud^] of 
A = 0.805 to the Burkert profile (|37|l . The asymptotic rotation velocity is taken to 
be Voo = c/1000. 

Numerical computations of the reduced deflection angle for the NTS halo were 
performed by assuming that the halo lens is half-way between the observer and 
the source, i.e. As/As = 1/2. The result is shown in Fig. |31where a and 8 are 
measured in arcsecond, but in general the units of 9 depends on the oberserver-lens 
distance (further details, see Refs. 1211 179J) . We recognize prominent features of the 
NTS lens due to the linear increasing mass (|17(l . There exists a limiting value of the 
reduced deflection angle of about a — 0.93 arcsecond. Hence, if the true angular 
position of the source (3 increases, nevertheless the halo lens deflects the light with a 
constant value. This behavior would change if the mass of the halo becomes finite. 

We can compare this result with the strong gravitational source of a BS for 
different potentials U, when the BS lens is again half-way between the observer 
and the sourcJ^Eill Observable differences are found depending on the choice 
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of the self- interaction. In the case of a simple mass term t/(|$|) = m 2 |$| 2 , the 
largest possible value of a is 23.03 degrees with an image at about 9 = n x 2.88 
arcsec where n — n(D Q \, ui) is the distance factor which is a function of the distance 
from the observer to the lens and the scalar field frequency the inverse of which 
is associated with the BS radius. For non-relativistic BS approximations, smaller 
angles will occur. The angle 9 of the image position can have very different orders 
of magnitude, depending on n. For example, n — 1 fixes 9 to be measured in arcsec. 
Under the assumption that the BS mass is 10 10 M Q , the distance D a \ is about 100 
pc. If the distance factor is n — 10~ 3 , then 9 is measured in milli-arcsec and the 
BS-lens is at about 100 kpc from the observer. 

In compar ison , halos modelled by the Burkert profile will not produce strong 
lensing as well^l Differences of our NTS model and the Burkert fit can be traced 
back to the fact that the metric corresponding to an asymptotically constant velocity 

|6o| 

t'oo in general exhibits a deficit angle similarly as a cosmic string . The linear mass 
density peaks at the location of the string, thereby producing a strong lensing. 

The properties of dark matter halos inferred from weak lensing""! provides a 
strong support for the existence of dark matter, wheras alternative theories of grav- 
ity such as MOND can almost be excluded. 
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